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A note on a piecewise-linear Dnffing-type system 

Gheorghe Tigan*, Alessandro Astolfi^ 


Abstract 

In [IT] it was suggested that the number of limit cycles in a piecewise-linear system could 
be closely related to the number of zones, that is the number of parts of the phase plane where 
the system is linear. In this note we construct an example of a class of perturbed piecewise 
systems with n zones such that the first variation of the displacement function is identically 
zero. Then we conjecture that the system has no limit cycles using the second variation of 
the displacement function expressed for continuous functions. This system can be seen as a 
feedback system in control theory. 
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1 Introduction 

In this work we consider a piecewise Dufhng-type system. The classical Duffing system deals with 
continuous functions while the system addressed in this paper is a discontinuous one. This system 
can be seen as a feedback system in control theory. Such non-smooth systems appear naturally in 
many practical systems because many physical phenomena presents discontinuities. For example 
in control theory, systems controlled by switching belong to this class of non-smooth systems. It is 
known that switching occurs in control systems in industry such as multi-body systems, intelligent 
systems, robots [T]. Neuronal firing in biology [5] or impacts in mechanics are other systems 
which address the problem. The first works on this topic are [3] and [4]. Recently results can be 
found in 0, i and [7]. In some cases, piecewise-linear systems offer a good approximations for 
nonlinear complex systems offering a valuable tool for investigating nonlinear phenomena. As it 
was pointed out in [7] , there is a feeling that piecewise-linear systems can present all the features 
met in nonlinear dynamics, such as homoclinic or heteroclinic orbits, limit cycles and attractors. 
In fact, Chua and collaborators discovered chaotic behavior in piecewise linear systems [5]. We 
will focus in the present work on the bifurcation phenomena related to the existence of invariant 
closed curves, such as limit cycles. The number and distribution (location) of limit cycles is one 
of the most important problems in qualitative theory of dynamical systems. There are numerous 
works on existence, number and distribution of limit cycles for continuous dynamical systems, for 
example [n m n m m [ig but not many exist for non-smooth systems. Investigations of the 
dynamics of non-smooth systems with one, two or three lines of discontinuity are performed in 
m, 0,0- Few papers deal with systems with a large number of discontinuities and not much 
is known about their dynamics. A recent work, which considers a piecewise-linear function for a 

*Department of Mathematics, ’’Politehnica” University of Timisoara, P-ta Victoriei, Nr.2, 300006, Timisoara, 
Timis, Romania, email: gheorghe.tigan@mat.upt.ro 

^Electrical and Electronic Engineering, Imperial College London, South Kensington Campus, London SWT 2AZ, 
UK, email: a.astolfi@imperial.ac.uk 


1 



Lienard system with n zones is m In that work it is suggested that the number of limit cycles is 
closely related to the number of zones. In the present work, we construct a piecewise Duffing-type 
system with n zones and conjecture it has no limit cycles. 


2 A piecewise-linear DufRng-type system 

Assume H is a Hamilton function and consider the following perturbed differential system 

x= dH/dy + ef{x,y,e) 

y= —dH/dx + eg{x,y,e) ^ ' 

where f{x,y,e) and g(x,y,e) are two enough smooth functions in x,y which depend analytically 
on a small parameter e. For any h on a real interval (a,6), we suppose that the set {(x,y) € 
M.;]I(x,y) = h} contains a closed curve Ch free of critical points (a circle for example) which 
depends continuously on h. Such a family of closed curves Ch corresponds to an annulus A of 
periodic solutions of the unperturbed Hamiltonian system dH = 0, that is, the system ([1]) for 
e = 0. 

Recall that, if we fix a transversal segment to the flow in o and parameterize it using the energy 
level h, then the function 

d{h,e) := P{h,e) - h = eMi{h) + e'^M 2 {h) + ... + e'^Mkih) + he {a,b), (2) 

where P{h,e) is the first return map (or Poincare map), is the displacement function defined for 
small e. The variations of the displacement function, Mk{h), are also called the Melnikov functions. 
Computing explicitly the Melnikov functions is a challenging problem and as long as we know they 
are determined only in some particular cases. It is known that the number of zeros of the first 
non-disappearing k — th order Melnikov function Mk(h) provides the upper bound of the number of 
limit cycles of the perturbed system emerging from the periodic orbits of the unperturbed system 
[18] . More exactly, if the first not identically null Melnikov function is Mk{h), then we have the 
following result: 

Theorem 2.1. If Mi{h) = ... = Mk-i{h) = 0, Mk{h) yf 0 for some h and hi is a root of the 
Melnikov function Mk{h) such that the m—derivative Mjf^\hi) ^ 0, m > 1, then for e ^ 0 
sufficiently small, system m has one limit cycle of multiplicity m in an 0(e) neighborhood ofCh^. 
In case that Mk(h) ^ 0 for any h, then for e ^ 0 sufficiently small, the system (Qp has no limit 
cycles in an 0(e) neighborhood of Ch. 

More details can be found in [2T| and m- The first Melnikov function for the system m , as it is 
reported in [22|, is given by: 

Mi(h) = (f g(x,y,0)dx - f(x,y,0)dy. (3) 

Jc{h) 

In [22] is presented a method for determining the second Melnikov function, following an al¬ 
gorithm described in [23] . for a system of type (|T|) with the Hamiltonian expressed in the form 
H(x,y) = ^y^ — U(x), where U(x) is a polynomial of degree at least 2, in the case when the first 
Melnikov function Mi(h) = 0. In this case, the second order Melnikov function is given by: 

M 2 (h) = ^(x,y,0)dx - ^(x,y,0)dy, (4) 
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provided that 


For similar systems of the form 


^(x,y, 0 ) + -^{x,y,0) = 0 . 


{ ^ \ \ 9i{x,y,e) \ 

\ y J V /V 92ix,y,e) ) ’ 


( 5 ) 


where {x,y) G 0 < e < 1 and / = (/i(x, y),/ 2 (x, ?/)), g = {gi(x,y,e), g 2 ix,y,e)) are two 
sufficiently smooth functions, the first Melnikov function is treated in [21] (used in [17]) and is 
given by: 


Mi{r)=[ e -^ 0 - 5 i(Tr(t), 0 )/ 2 (Tr(t)))dt ( 6 ) 

Jo 

where we assumed that the system ([5]) for e = 0 possesses a family of periodic orbits Cr : Tr{t),r > 0 
depending on a positive real parameter r and having the period Tr. This formula can be used even 
though / is differentiable but g is discontinuous in some isolated points, as remarked in [17]. It 
will be employed in the present work in such a case. 

Consider now the planar discontinuous Duffing-type system given by: 


X 

y 


-1 0 ) (j )■"'( 9 ^ 9 .^) ) 


( 7 ) 


where 5 is a real function defined on the set J x R x (—ei, £ 1 ), with I = (—00, oi] U (ai, 02 ] U ... U 
(a„, 00), Oi G R, f = 1, 2,..., n, £1 > 0, discontinuous in the points at and, for a fixed £, differentiable 
with respect to (x,y) on each of the (n + 1) strips (ai,ai+i) where it is defined. If g{x,y,e) = x^ 
and e = f3 we meet the classical Duffing system [H] in a particular case. In [^ we studied a 
Duffing continuous system identifying the conditions to transition to chaos. 

Our first result from this paper is stated in the following Proposition: 


Proposition 2.1. Consider two increasing sequences of real numbers 

ao = —00 < 0 < tti < ... < a„ < a„+i = +00, —00 < ao < cti < ... < an < +00, n G N, n > 1 and 
the non-smooth linear function with n + 1 zones 


g{x,y,s) = aiX + ey */x G (oi, Oi+i], f = 0,1,..., n ( 8 ) 

with the convention that the last interval is (ori,a„+i). 

Then, for any 0 < £ ^ 1 and any two sequences (ai)i=o,i,...,n+i and (ai)i=o,i,...,rt as above, the 
first Melnikov function of the system & is identically zero. 

PROOF It is clear that the system © is of type ([5]) with /i = y,f2 = -x,gi = 0 and 
g 2 {x,y,e) = g{x,y,e) = a^x + ey. With these functions, the unperturbed system ([7]) has a one- 
parameter family of periodic solutions which are circles of the form Cr '. x^ -\-y^ = r^. Choosing the 
parametrization Tr{t) : x = rsint, y = rcost, t G [0, 27 r), and denoting g{x) := g 2 (x,y,Q) = a^x, 
then using m we get 


Mi(r) 


r2Tr 


rcost g{r smt)dt. 


(9) 
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One can observe that the same result is recovered from the formula m, at least formally since g 
is discontinuous, applied for systems of the form o with H{x,y) = f{x,y,e) = 0 and 

g{x,y,e) = UiX + ey, denoting the circle Cr ■ x'^ + y^ = and g{x) := g{x,y,0) = aix. Indeed, 
from we have 

Mi{r) = j) g{x)dx, (10) 

and using the same parametrization x = r sint, y = r cost, t £ [0, 27r), it leads to ([51)• 

Find now Mi{r). From ([21) we have 

p27T 

Mi(r) = / r cost g{r sin t)dt (11) 

^0 

/2 ^37r/2 p2TT 

r cost g{r s\nt)dt + / r cost g{r smt)dt + / r cost g{r sint)dt. 

Jtt/2 J3it/2 

Compute in the following the three integrals. Let m £ {0,1, ...,n} such that am < r < 0^+1 
and to = 0 < ti < t 2 < ■■■ < tm < tm+i = 7’‘/2 an increasing sequence of real numbers in 
[0,7r/2] given by sint^ = 2±^ i = 1,2, It is obvious now that if t £ (ti,ti+i), then rsint £ 

(oi, Oi+i), i = 0,1, 2,..., TO. Consequently 



r'iT/2 


m 

'ti 


M^{r):= / r cost g{r sin t)dt = 'y^ / cost sin tdt 

•^0 i =0 

J.2 ^ 

= Y X! - sin^ U) 

^ m 

- -ai-i)al 


2=0 

.2 1 


= 2 


2=1 


( 12 ) 

(13) 


For the second integral, consider a decreasing sequence of real numbers tm+i = 7r/2 < tm < 
tm-i < < ti < to = 37r/2 given by sint^ = —, t = 1,2 ,...,to. Because the sin function is 

decreasing on the interval [7r/2,37r/2], one gets that, if t £ (ti+i,ti), then rsint £ (ai,ai+i),t = 
0,1, 2,..., TO. Therefore 


|•3^r/2 rn 

M^{r):= / r cost 5(r sin t)dt = ^ cost sin tdt (14) 

JTr/2 Jti+1 

J.2 

= y ^ai(sin^ti - sin^t^+i) (15) 

i=o 

r2 1 ™ 

~ y (®o Oi-m) ~b 2 ^ ^ (*^» ~ 0!i—i)a^. 

^ i=l 

In the last case, because rsint < 0 for t £ [37r/2, 27r) we have that rsint £ (ao,ai) so that 

5(rsint) = oorsint. 

Therefore 
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rcost g{rsmt)dt = 



r^ao costsintdt = ——ao. 


(16) 


Finally, 


Mi(r) = M^{r) + M‘^{r) + M^{r) = 0. 


(17) 


In the following we consider a much more general class of discontinuous functions 

g{x,y,e) = aih'{x) + ey ifx£{ai,ai+i],i = 0,l,...,n (18) 

with the same convention for the last interval, i.e. it is (a„, o„+i), where h'{x) is the derivative of 
a differentiable function h(x) and we will prove a similar result given by: 

Proposition 2.2. If g is a non-smooth function given by m , then, for any 0 < e <C 1 and any 
two sequences {ai)i-o i „+i and (ai)i=o i n as above, the first Melnikov function of the system 
& is identically zero. 

The proof is similar. Proceeding as above, we have that: 

^7r/2 rn 

M^{r):= / r cost (/(r sin t)dt = ^ / rai cost h'{r smt)dt (19) 

m 

= h(r)am - /i(0)ao - “ ai-i)h{ai), 

i=l 

/•7r/2 m .ti 

M‘^{r):= / r cost g{r smt)dt = '^ / rai cost h'{r sin t)dt (20) 

•to Jti+i 

m 

= -h{r)am + h{-r)a^ + X]*'*^* “ ai-i)h{ai), 

2=1 


and 

p2-K p2-K 

M^{r) •.= / r cost g{r s\nt)dt = / rao cos th'{r sin t)dt = h{{))aQ — h{—r)aQ, (21) 

J'iTT/2 Js-k/2 

so we arrive to the same result: 


Mi{r) = M^{r) + M'^fr) + M^{r) = 0. 


( 22 ) 
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As the first Melnikov function is identically zero, we can say nothing about the number of limit 
cycles. However, we conjecture the following fact: 

CONJECTURE: The non-smooth system & with the function g given by m has no limit cyeles. 

We base this conjecture of the following fact. Compute the second Melnikov function for the 
system Q with g given by (fTSll . One can observe that ^{x,y,0) + ^{x,y,0) = 0 on any strip 
where g is defined, so using o we get that: 

M 2 (r) = </ ^{x,y,0)dx — ^{x,y,0)dy = ydx = [ cos^ tdt = ^ 0. (23) 

Jcr Jcr do 

As M 2 (r) = 0, r > 0 has no root the above Conjecture is justified by Theorem 12. II However, we 
can not present this Conjecture as a result, because the formula of the second Melnikov function 
which we used from [22| was proved for continuous functions but we applied it for functions with 
a finite number of discontinuities. 


3 Conclusions 

In this work we started to explore an example of a class of non-smooth dynamical systems with the 
first Melnikov function identically zero. We applied the formula of the second Melnikov function 
reported in a work of Iliev [22j in order to investigate further the existence of limit cycles. This 
formula has been deduced for continuous functions but it is quite plausible it remains valid for 
piecewise functions. 
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